Hole concentration in a diluted ferromagnetic semiconductor 
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We consider a mean-field approach to the hole-mediated ferromagnetism in III-V Mn-based semicon- 
ductor compounds to discuss the dependence of the hole density on that of Mn sites in Gai-^Mn^As. 
The hole concentration, p, as a function of the fraction of Mn sites, x, is parametrized in terms of the 
product rn* Jp d (where m* is the hole effective mass and J p d is the Kondo-like hole/local-moment 
coupling), and the critical temperature T c . By using experimental data for these quantities, we 
have established the dependence of the hole concentration with x, which can be associated with the 
occurrence of a reentrant metal-insulator transition taking place in the hole gas. We also calculate 
the dependence of the Mn magnetization with x, for different temperatures (T), and found that as 
T increases, the width of the composition-dependent magnetization decreases drammatically, and 
that the magnetization maxima also decrease, indicating the need for quality-control of Mn-doping 
composition in diluted magnetic semiconductor devices. 



PACS Nos. 71.55.Eq, 75.30.Hx, 75.20.Hr, 61.72. Vv, 72.25.Dc 



Over the last few decades, a considerable amount of 
work has been devoted to the understanding of elec- 
tronic, optical and transport properties of diluted mag- 
netic semiconductors (DMS). Interest in these materials 
was boosted in the early 1990's with the discovery of 
ferromagnetism in III-V materials alloyed with transi- 
tion elements like Mn |l|,||. Ferromagnetic semiconduc- 
tors bring about the possibility of controlling both spin 
and charge degrees of freedom, which, when combined 
with the capability of growing low-dimensional struc- 
tures, opens up exciting new prospects for the production 
of spintronic devices. Potential applications include non- 
volatile memory systems and quantum computing 

Especial attention has been focused on Gai-^Mn^As 
alloys, which exhibit very interesting magnetic and trans- 
port properties. Mn atoms have five electrons in the 3d 
levels and two electrons in the 4s levels, and their incorpo- 
ration into a GaAs matrix plays two roles: they act both 
as S = 5/2 local moments, and as acceptors generating 
hole states in the material. The equilibrium solubility of 
Mn atoms in GaAs is quite low, being only of the order of 
10 19 cm -3 j|. However, with the use of molecular-beam 
epitaxy techniques at low temperatures, several groups 
have recently succeeded in producing homogeneous sam- 
ples of Gai-^Mn^As with x as high as 0.071. It has been 
observed that for 0.015 < x < 0.071 the systems become 
ferromagnetic, with doping-dependent critical tempera- 
tures T c (x) reaching a maximum of 110 K for x — 0.053 

The appearance of a ferromagnetic state in these ma- 
terials has been attributed to an exchange coupling be- 
tween the localized Mn moments mediated by the holes, 
whose strength should depend on the hole concentration 
p. In principle, one would expect that each Mn would 
provide one hole, leading to a density of holes equal to 
that of the magnetic ions. However, while an accurate de- 
termination of the hole concentration is hindered by the 



anomalous Hall term, experimental data indicate that p is 
only a 15 to 30% fraction of that of magnetic ions [p~0| p~3|| . 
The mechanism responsible for the discrepancy between 
hole and Mn densities is not clear. As pointed out by 
Matsukura et al. JTlJ , such discrepancy might be due to 
compensation of Mn acceptors by deep donors such as 
As antisites, which are known to be present at high con- 
centration in low-temperature grown GaAs JOj . Another 
possibility would be the formation of sixfold-coordinated 
centres with As (Mn 6As ), which would compensate Mn 
atoms on substitutional Ga lattice sites As a con- 

sequence, the relation between hole concentration and 
that of Mn has not been so far theoretically established, 
which would be of great interest for the design of new 
devices. Our main purpose here is to present a quantita- 
tive analysis on this issue, based on a simple model for 
the magnetic behaviour of these systems. 

We adopt the generally accepted view that a given 
Mn ion interacts with the holes via a local antiferromag- 
netic Kondo-like exchange coupling J pc [ between their 
moments ll^Jl^ -|l9|] . This interaction is thought to lead 
to the polarization of the hole subsystem, which would 
then give rise to an effective ferromagnetic coupling be- 
tween the Mn moments. Though there has been some 
debate as far as the details of the above picture are con- 
cerned (e.g., whether or not such effective interaction is 
well described by an RKKY term [gO 21 1), there is an 
overall consensus on the fundamental role played by the 
hole-mediated mechanism. At any rate, the approach we 
follow here does not depend on the details of the effective 
Mn-Mn interaction. 

We start with a Hamiltonian for the two coupled sub- 
systems of the form 

n = Hun + n h + j pd s i • s * 6 ( r * - R /) - (!) 

i,I 

where 7iMn describes the direct {i.e., non-hole- mediated) 
antiferromagnetic exchange between Mn spins, 7ih de- 
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scribes the hole subsystem, and the last term corre- 
sponds to the aforementioned Mn-hole exchange inter- 
action, with Si and s, labeling the localized Mn spins 
(S = 5/2) and the hole spins (s = 1/2 ), respectively. As 
a first approach, we neglect Ttyin and consider Hh within 
a parabolic-band effective-mass approximation; we com- 
ment below on more general descriptions of TL^. 

Within a mean-field approximation, the Mn magneti- 
zation is given by 



M = N Mn gii B Mi = nungmSBs 



JpdS 

2k B T 



M h 



(2) 



where riMn = xn s is the density of Mn ions, with n s being 
the density of Ga lattice sites, Mi is the magnetization 
density of the Mn ions, g = 2 is the Mn Lande g-factor, 
and E>s[- ■ .] is the Brillouin function. The magnetization 
density of the hole subsystem, Mh = (n-[ — "j,), is sup- 
posed to be uniform within the length scale of magnetic 
interactions, so it can be calculated self-consistently by 
considering a Fermi sea of holes with effective mass m* , 
in the presence of the mean magnetic field generated by 
the Mn ions; it is therefore given by 



Mh = A J pd x Mi p 

rrip 



where A = 6 (1/3tt 2 ) 2/3 (m e /h 2 a 3 ) , 
electron mass; a = 5.65 A is the GaAs lattice constant. 

The critical temperature as a function of the hole den- 
sity and of the Mn composition is obtained by linearizing 
the self-consistency relations given by Eqs. (0) and (0): 



1/3 



(3) 



is the free- 
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S(S + 1) [{m* /m e )J 2 pd ] xp 1 ' 3 



(4) 



Specializing Eq. (Q) to S = 5/2, we write the hole con- 
centration as 



T c (x) 



(m*/m e )J 2 d 



(5) 



where £ = 5.29 x 10 , in units such that J pd is given 
in eV nm 3 . 

In view of the uncertainty on the available experimen- 
tal values for m* and J p d, and to the difficulties in obtain- 
ing accurate estimates over a wide range of hole densities, 
p, the following strategy is adopted. We first use the fact 
that Hall resistance measurements p2| yield an unam- 
biguous |2j| value of p — 3.5 x 10 20 cm -3 for the sample 
with x = 0.053, for which T c = 110 K. We then take 
these values into Eq. (||), to fit a value for the product 
(m*/m e )J 2 d = 2.4 x 10~ 2 (eV nm 3 ) 2 . And, finally, we 
use this value, together with the experimental transport 
data p0[ | for T c (x), to obtain p over a wide range of x, 
shown as filled squares in Fig. |[ The error bars in Fig. |l| 
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FIG. 1. Theoretical results for the hole concentration as a 
function of the fraction of Mn sites for a (Ga,Mn)As ferro- 
magnetic alloy. The dashed line corresponds to a hole con- 
centration equal to that of Mn sites, whereas the filled squares 
are the present mean-field result; the full curve is obtained as 
reasoned in the text. I and M respectively denote insulating 
and metallic phases. 



reflect the uncertainties in the determination of T c (x), as 
displayed in Fig. 3(c) of Ref. 113]. The adequacy of this 
procedure is illustrated in FigT[lJ The calculated values 
for p(x) lie below the concentration of Mn ions, shown 
as a dashed line, in agreement with experiment. We also 
highlight in Fig. [j] the boundaries of the metal-insulator 
transitions (MIT's), as determined from resistivity mea- 
surements |10| . The present theoretical estimates for p in 
the insulating phases are based on the assumption that 
the localization length in insulating samples, though fi- 
nite, is significantl y la rger than the length scale of mag- 
netic interactions |l3[ , in which case the present mean- 
field approach is a good starting point. 

Before accepting these estimates for p{x) at face value, 
one should note that a closer look at the experimental 
data for T c (x) suggests a linear behaviour in the 
range of x of the order 0.015-0.035 which would imply, 
through Eq. (JsJ) , a constant p in that range; this constant 
behaviour, however should not prevail at low concentra- 
tions, x — > 0, and presumably one should havep oc x — > 0. 
These considerations are incorporated in the full curve 
displayed in Fig. [l], which lies within the error bars of 
the calculated hole concentrations. 

Our theoretical estimates for p(x) are therefore 
strongly suggestive of p(x) reaching a maximum value 
within the metallic phase. As a consequence, all attempts 
to increase T c should be carried out for samples in the 
metallic phase, for Mn concentrations about 0.05. More- 
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FIG. 2. Schematic density of states (DOS) versus energy 
for the impurity band, for up- (top) and down-spins (bot- 
tom). Under each DOS curve, the hashed and empty regions 
correspond, respectively, to delocalized and localized states; 
these are separated by mobility edges. The exchange splitting 
is proportional to the off-set between the j and J, bands, and 
the Fermi energy (ef) increases to the right, towards the top 
of the valence band (not shown). 

over, notwithstanding the considerable uncertainties p^ j 
in the measurements of p, the data shown in Fig. [I] are 
in qualitative agreement with those obtained from Hall 
measurements by Matsukura et al. JO ]; as we discuss 
below, this is also consistent with findings from recent 
photoemission spectroscopy measurements p5| . 

At this point several comments are in order. First, the 
model is indeed very simple, for it does not incorporate 
aspects such as a Kohn-Luttinger treatment of the va- 
lence states [|13|J26| , effects of impurity potentials, a site 
energy term arising from the Mn potential, a correlation 
energy representing hole-hole repulsion, and so forth. In 
addition, the model is treated within a mean-field ap- 
proximation which neglects fluctuations in spin, charge, 
and disorder degrees of freedom. Nonetheless, one ex- 
pects these limitations to be minimized, to some extent, 
by the fact that experimental data for T c (x) are used as 
input. However, it is exactly this simplicity that allows 
us to obtain a direct relation between hole concentration 
and Mn fraction, which, in turn, can be promptly used 
as a rough guide to experiments. Clearly, the present 
results must be viewed as a first approximation to p{x), 
since one should still be able to obtain such relation phe- 
nomenologically through improved models and approx- 
imations, though with a considerable amount of extra 
computational effort. Consider, for instance, the case of 
Monte Carlo simulations of the Kohn-Luttinger Hamilto- 
nian for the semiconductor valence bands p7| : since the 
dependence of T c on J pc i and on m* is different from that 
of the mean-field prediction [Eq. (Q)], the one-parameter 
fitting strategy adopted here is not so straightfowardly 
applicable; these parameters would have to be separately 
adjusted, demanding many additional runs. It is there- 
fore hardly surprising that most of the improved theo- 
retical approaches |l3| , |27| -p9| consider p (instead of x) as 
an independent variable and, accordingly, present plots 
of T c (p), for fixed x; the issue of the relation between p 
and x is then set aside. 

The present approach also leads to a qualitative under- 
standing for the dependence of p with x being essentially 
related to the occurrence of MIT's taking place in the hole 
subsystem. Within our approximation, the Fermi energy 



tracks the behaviour of p, since £f oc p 2 ^ 3 , while the ex- 
change splitting A oc x. Figure ||(a) shows the schematic 
impurity bands for each spin channel, in the very-low 
doping regime in which the gas is supposed to be unpolar- 
ized. As x increases, the gas can sustain polarization and 
still be insulating, provided the Fermi energy lies below 
the mobility edge, as shown in Fig. |^(b). Further increase 
in x causes £f to increase and to lie within the delocal- 
ized states of the up-spin impurity band, as depicted in 
Fig. ||(c): The system becomes metallic. Whether or not 
the Fermi energy also lies within the delocalized region of 
down-spin impurity band is a very interesting question, 
which cannot be answered by our simple model; the solu- 
tion of this particular issue should have bearings on the 
efficiency of Gai-^Mn-j As-based devices as spin filters. 
Once £f reaches a maximum within the metallic phase, 
its initial decrease upon increasing x is compensated by 
an increase in A, so that the Fermi level still lies within 
the delocalized states. However, with continuing increase 
in x the exchange splitting can no longer make up for 
the decrease in £p, and the latter eventually crosses the 
mobility edge again, lying within localized states [Fig. 
0(d)] : The system reenters an insulating phase. One 
may argue that a description in terms of impurity lev- 
els rather than impurity bands may be more adequate in 
the range of Mn concentrations considered here. Even so, 
the movement of the Fermi energy described above is still 
applicable with slight modifications: the metallic phase 
would then correspond to Sf reaching the top of the va- 
lence band. This latter picture is actually in qualitative 
agreement with recent photoemission measurements pq] 
of the Fermi level as a function of Mn concentration in 
Mn^Gai^^As. 

We now discuss the magnetization of the Mn ions, as 
obtained by solving Eqs. (g) and (||) for M for a given 
Mn composition and temperature. The mean-field the- 
oretical results are shown in Fig. || for various temper- 
atures; for T = 75 K and 100 K, we had to resort to 
an interpolation of the experimental T c data from Ohno 
and Matsukura JTcp ■ Two effects are apparent from the 
calculated results. Firstly, the magnetization maxima 
with respect to x decrease with increasing temperature, 
as it would be expected since one approaches the crit- 
ical temperature from below [see also Fig. 3(c) of Ref. 
fllPf]. Moreover, the widths of the composition-dependent 
magnetization curves decrease quite drammatically with 
temperature. An immediate consequence of these results 
is that DMS (Ga,Mn)As device applications at tempera- 
tures < T c would require a definite quality-control of the 
Mn doping composition. 

In summary, we have established a theoretical sce- 
nario for the behaviour of the hole concentration in 
Mn^Gai^As as a function of both x and T c , based on 
a simple mean-field approximation to the hole-mediated 
ferromagnetic Hamiltonian. In our picture, the concen- 
tration of holes is approximately constant in the low- 
doping insulating phase, then rises to a maximum in the 
metallic phase, and drops again in the reentrant insu- 
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FIG. 3. Mn contribution to the magnetization [Eq. (^|)] as 
a function of the fraction of Mn sites, for four different tem- 
peratures, T = 40 K, 50 K, 75 K, and 100 K. 



lating phase. Our approach also allows one to view the 
underlying mechanism of the reentrant MIT's as an oscil- 
lation of the Fermi energy caused by a delicate balance 
between band filling and exchange-splitting. We have 
also noted that the larger the temperature, the narrower 
is the range of compositions leading to a non-zero Mn 
magnetization. The present approach should certainly be 
extended to include a more complete description of the 
acceptor states, taking into account the spin degrees of 
freedom, spin-orbit coupling, compressive/tensile strains, 
etc. Moreover, a proper treatment of disorder - e.g., by 
explicitly considering a random, instead of continuous, 
distribution of Mn ions - should lead to a more realistic 
description of the MIT. In this respect, many-body effects 
due to correlation among the holes should also influence 
p(x), especially in the insulating phase. Of course, an ap- 
propriate description of the physical mechanisms related 
to As antisites and Mn 6As centres, as far as the hole vs. 
Mn concentrations is concerned, is certainly a formidable 
task, which nonetheless deserves future theoretical atten- 
tion. 
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